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Abstract 

In this paper we calculate, in the large N limit, the eigenvalue 
density of an infinite product of random unitary matrices, each of 
them generated by a random hermitian matrix. This is equivalent 
to solving unitary diffusion generated by a hamiltonian random in 
time. We find that the result is universal and depends only on the 
second moment of the generator of the stochastic evolution. We find 
indications of critical behavior (eigenvalue spacing scaling like 1/iV 3 / 4 ) 
close to 9 = 7T for a specific critical evolution time t c . 



1 Introduction 

A key feature of random matrix theory is that many properties of random 
matrix models do not depend on the fine details of these models but only 
on some very general symmetry properties and a very limited number of 
numerical coefficients (usually just a single coefficient is enough). These 
universal properties facilitate the widespread applications of random matrix 

*e-mail: uf r j anikOif . uj . edu . pi 
t e-mail: wieczor@th.if.uj.edu.pl 



1 



models in various fields since one can use the models to learn something about 
the behavior of complex systems without knowing all the precise microscopic 
details of these systems. 

The simplest example of such a behavior is the eigenvalue density of a 
matrix with entries independently distributed according to some probabil- 
ity distribution. Then the eigenvalue density in the N — > oo limit follows 
Wigner's semicircle law with the scale set just by the second moment of the 
distribution. All the dependence on other properties of the initial probabil- 
ity distribution disappears. In this work we will uncover a similar type of 
universality in a different context. 

The aim of this paper is study infinite products of unitary random matri- 
ces, and derive their properties in the large N limit. This can be interpreted 
as a multiplicative diffusion process on the unitary group 1 . A natural physi- 
cal interpretation would be of a quantum-mechanical evolution governed by 
a hamiltonian which changes randomly in time. Another possibility would 
be the modelling of Wilson loops in lattice gauge theory. In this paper we 
will not examine further these possible application but rather concentrate on 
mathematically solving the model. 

We show that the eigenvalue density exhibits universality properties i.e. 
it only depends on the second moment of the random hamiltonian which 
generates the stochastic evolution. But of course the resulting eigenvalue 
density is much more complex than the semicircle law. We derive equations 
for the eigenvalue density, give explicit expression for the lowest moments 
and study the properties of the model close to a phase transition where the 
eigenvalue support begins to cover the whole unit circle. 



2 Multiplicative unitary diffusion 

We consider the product of M unitary N x N random matrices Uk in the 
M, N -> oo limit: 

M 

U = lim lim FT U k , (1) 

where the C4's are generated by 

U k = e UH \ (2) 

1 Recently matrix valued multiplicative diffusion has been considered for 2 x 2 real 
matrices in [1] and for infinite hermitian and complex matrices in [2]. 
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and where e = \Jt/M. Such a scaling is standard for diffusive processes and 
works also very well for matrix-valued diffusion processes studied in [1] and 
[2]. t is then a real parameter corresponding to 'diffusive' evolution time and 
the continuum limit M — > oo exists. The generators of the evolution are 
N x N hermitian matrices drawn from a probability distribution 

P(H) ~ e - Ntrv ^ (3) 

where we assume that the first moment mi = ^ -^tr if ^ vanishes mi = 0. We 
will show below that the spectral properties of (1) depend only on the second 
moment m 2 of the distribution (3) 

m 2 = (±-trH 2 ). (4) 



,N 

The main aim of this paper is to find the eigenvalue distribution of the 
product (f) 



f N 




ti0,t) = l-y:ie«8(e«-e«>)) , (5) 



where the e %0j are the eigenvalues of U defined through (1). 

In the next section we will use free random variable methods to derive an 
equation from which one can get p(9,t). Sometimes we will omit the second 
argument but of course the dependence on t will be there. 



3 The S-transform method 

The main difficulty encountered in [2] when considering products of random 
matrices was the necessity to deal with nonhermitian matrices and eigenval- 
ues covering two dimensional regions of the complex plane. Here fortunately, 
since the product matrices are always unitary, the eigenvalues lie on the unit 
circle and hence they can be uniquely reconstructed from just the knowl- 
edge of the moments. Therefore all information is encoded in the asymptotic 
expansion of the Green's function 

Gw=rJ%». (6) 

Jo z — e w 

Our aim now is to obtain the spectral density Pprod{0) of product U = 
YlkLiUk, equivalently the corresponding Green's function Gprod(z), from 
the spectral density Ph(Q) f° r the generator H. 
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To do that we use S-transforms introduced in [3]. Firstly we define an 
auxiliary function x( z ) through: 

-g(-)-1 = z; (7) 

X \XJ 



Then the S-transform is: 



S(z) = — x(^) ■ (8) 



Their main property is that the S'-transform of a product of random matrices 
is a product of S'-transforms of the individual factors. Before we apply this 
setup to (1) let us note that putting together the two previous equations we 
arrive at a functional relation satisfied by S and G: 

In our case all single matrix Green's functions are the same since they 
come from the same distribution so we can write: 

M 

Sprod(z) = lim I] Si(z) = hm (S 1 (z)) M , (10) 

M— >oo. , M—KZD 

1=1 

Let us now find S 1 (z). The eigenvalue density of Lq = e ieH is given by 

Pi{e,e) = \ PH (°-\ . (11) 



where e = yt/M. Inserting (11) and the definition (6) into (9) leads to an 
equation for Si: 

- - (') 5 M = 1. (12) 

From the form of (10) we see that we need to calculate Si only to the order 
0(e 2 ). Substituting u = 9/e and expanding in Taylor series in e we obtain: 

j p H {u) (l + ieuz + (-wV - ^u 2 z + sz)e 2 + 0(e 3 )^j du = 1 , (13) 

where s = s(z) comes from Taylor expansion 2 Si(z) = 1 + s(z)e 2 + 0(e 3 ). 
2 Here we used the assumption that the first moment of H vanishes. 



4 



We may now calculate s(z) 



s(z) = (z + 1/2) (u 2 ) = (z + l/2)m 2 . (14) 

From (10) we may now obtain the S'-transform for product: 

/ t \ M 

S PROD = lim (S 1 ) M = lim 1 + —s(z)) = e« z+ v*» . (15) 

This result shows that Sppon depends only on the second moment m 2 
of H. So in the limit M — > oo we obtain universal behavior of the system 
independent of the spectral density of generator of the stochastic evolution 
H, as long as the first moment vanishes (no drift) and the second moment 
is finite. It will be interesting to consider cases where these assumptions are 
violated which would lead to anomalous diffusion. We leave these problems 
for future investigation. 

The final step is to come back from SppoD to the Green's function 
GpROD(z,t) (from now on we drop the subscript). It is convenient to in- 
troduce auxiliary function f(z,t) as: 

G(z,t) = 1 -±I^. (16) 



It is easy to check that / fulfills an equation: / (^yj — z. This means 
that / and 1/x are functional inverses of each other so the following relation 
Vx(/) = z is also true. This observation, together with the result (15) and 
the definition (8) leads us to the final equation: 

zf = {1 + f) e-W+D"* . (17) 

This equation encodes all the spectral properties of the unitary diffusion 
process (1). In the next section we will proceed to investigate some of its 
properties. 



4 The dynamical properties of the unitary 
diffusion 

In this section we analyze the dynamical behavior of unitary matrix diffusion. 
For small times t the eigenvalues will be concentrated only in a small neigh- 
bourhood of 6 = 0. For longer times the support of the eigenvalue density p 
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Figure 1: The time evolution of the spectral function p(9,t) (the dots 
represent numerical simulation). The figures show p(9,t) after a time 
t = 1,250,500 up to 2750 (with m 2 = 1/500). The solid lines represent 
the eigenvalue density (29) obtained by the numerical solution of (17). 



will expand and when some critical time t c is reached the eigenvalues will fill 
the whole circle. But of course the eigenvalue density p will be nonuniform. 
In fact we expect a critical behavior close to 6 = n with nonstandard frac- 
tional eigenvalue spacing. Only later for t — > oo the eigenvalues will become 
uniformly spread over the whole unit circle. In fig. 1 we show numerical 
results for the eigenvalue density obtained by generating unitary matrices 
and compare it to the one extracted from (17) (see below). 
In this section we will quantitatively analyze this behavior. 
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The support of the eigenvalue distribution and the crit- 
ical time t c 

Although one cannot find an analytical formula for the eigenvalue density one 
can analytically find the edges of the eigenvalue support. These occur when 
the Green's function has an infinite derivative d z G = oo. Differentiating (17) 
with respect to z gives 



1 = -j[ (l + m 2 tf + m 2 tf) e -*(/+D™2 . (i 8 ) 



So the end points are determined through the solutions of equation l+m 2 tf+ 
m 2 tf 2 = 0. Once we know / we can reconstruct the end-points z using (17). 
The result is 

_ yg - m 2 t + i^t ^ tV ^- t 
Zed9e V^^i-t^i' 6 1 yJ 

and its complex conjugate z* edge . When these two solution are equal (z edge = 
z tdge = ~~ l)^ t ne eigenvalues will cover the whole circle. This will happen for 
the critical time 

t c = — ■ (20) 

The moments of U 

Another quantity which can be analytically calculated are the moments of 
U . The coefficients of the auxiliary function / around z = 00 



00 



/(*) = E^> (21) 

k=i z 

are indeed directly linked to the moments: 

a k = ^trU k y (22) 

So inserting (21) into (17) allows us to find the moments iteratively. The 
expressions for the lowest ones are: 

cii = e~^ m2t (23) 
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a 2 = e- m2i (-l+m 2 t) (24) 

1 3 mo t r» o 

a 3 = - e — ~ (2-6m 2 t + 3m 2 t 2 ) (25) 
a 4 = -ie" 2m2 *(-3 + 18m 2 t-24m2t 2 + 8m^ 3 ) (26) 

] 5 mo t on on A A 

a 5 = — e — ^(24- 240m 2 t + 600m 2 t 2 - 500m 3 t 3 + 125m^ 4 )(27) 

In fig. 2 we compare the formulae for the lowest 4 moments with numerical 
simulations of the unitary matrix diffusion and find complete agreement. 

The eigenvalue density 

The equation (17) allows us to directly reconstruct the Green's function. 
However it is very simple to recover also the eigenvalue density. This follows 
from the observation that the moments of U are just the Fourier coefficients 
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of the eigenvalue density p(9) 

(^ Uk )= £ K K0)e M . (28) 

Using the relation of / to moments derived earlier, and the symmetry p(9) = 
p(—0) one finds finally 

?m = -h Re (l +f )- (29) 

Equation (17) may be easily solved numerically. In fig. 1, we show the 
resulting eigenvalue density together with numerical simulations for various 
times t. 



Critical behavior at t = t c and level spacing 

At t — t c the edges of the eigenvalue support touch at z — — 1. Typically in 
such cases we expect new critical type of behavior and nonstandard scaling 
of eigenvalue spacing with N. Let us analyze now this behavior. Inserting 
t = t c and / = -1/2 + F to (17) we obtain: 

To find the behavior close to z — — 1 (equivalent to 6 = n) we expand the 
left hand side of (30) in F and put z = — 1 + iy to get 

16F 3 

-l + iy (3i) 

Using the relation between / and the eigenvalue density (29) we thus find 
the behavior close to 9 = n: 



(32) 



Such behavior of the eigenvalue density leads to nonstandard eigenvalue 
spacing and signifies the appearance of new universal regime on the scale 
of eigenvalue spacing (analogous to Airy universality and l/A 2 / 3 spacing 
on the edges of the eigenvalue distribution of a generic hermitian random 









(l6y 


> %os i r 


|0-7T| 









9 



matrix [4, 5, 6, 7] in contrast to the standard 1/JV spacing in the classical 
Wigner-Dyson regime [8]). 

In our case, the number of eigenvalues between n and A is approximately 
equal to n ~ iV(A — 7r) 4//3 . Reexpressing A in terms of n shows that the 
eigenvalue spacing in the vicinity of 9 = n scales like A similar 

scaling appeared in a certain class of chiral random matrix models at finite 
temperature [9]. It would be interesting to compare these regimes and/or try 
to apply the methods of [10] to the case at hand. We leave this problem for 
further investigation. 

5 Discussion 

In this paper we considered multiplicative unitary matrix diffusion generated 
by random hermitian matrices. We found the eigenvalue density as a function 
of evolution time in the large N limit using S'-transform methods. The 
eigenvalue distribution turns out to be universal and depends only on the 
second moment of the random hermitian matrix which generates the diffusion 
process. 

We found that at a critical time of evolution t = t c the eigenvalues start 
to fill the whole unit circle, and close to = tt a nonstandard eigenvalue 
spacing ~ sets in which signifies the appearance of a new critical 

regime. 

There are various further issues that one could investigate. Firstly, re- 
laxing the assumption of the existence of the second moment might lead to 
defining anomalous diffusion processes. Secondly it would be interesting to 
study microscopic properties of these unitary matrices, however in order to do 
that new methods have to be developed. Thirdly a more detailed investiga- 
tion of the critical behavior at t — t c close to 6 = n would be interesting and 
last but not least the application of these results to some physical situations. 
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